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Abstract. Timescales for mixing of midlatitude air into the tropical lower 
stratosphere are deduced from observations of long-lived tracers N 2 O and CCI 3 F. 
Bulk mixing between tropical and midlatitude regions is assumed to be isentropic 
and relatively slow compared with local mixing within each region. The mean 
value of the mixing timescale ranges from 12 to 18 months near 20 km. There is a 
tendency for shorter mixing times at higher and lower altitudes, although vertical 
profiles of mixing cannot be definitively established by the data. A more robust 
quantity is given by the fraction of midlatitude air entrained into the tropical 
upwelling region. Implied mixing fractions exceed 50% above 22 km. 


Introduction 

It has long been recognized that the observed dis- 
tributions of stratospheric water vapor and ozone im- 
ply upward transport across the tropical tropopause 
and a mean meridional circulation characterized by as- 
cent in the tropics and subsidence at higher latitudes 
[Brewer, 1949; Dobson, 1956]. Early work by Newell 
[1963] on the dispersal of 185 W bomb debris and more 
recent studies of sulfate aerosol distributions [Trepte 
and Hitchman, 1992] show the persistence of maxima in 
tracer fields within 20° of the equator and between 21 
and 28 km altitude. These observations indicate a bar- 
rier to mixing with the extratropics, at least in compari- 
son to rapid irreversible mixing within midlatitude “surf 
zones” [McIntyre and Palmer, 1984]. Poleward trans- 
port across the subtropical barrier appears to be associ- 
ated with planetary wave activity, such as the tongues 
of air stretching from the tropics into midlatitudes re- 
ported by Randel et al. [1993]. These features have been 
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identified also in transport model simulations [Waugh, 
1993; Chen et al } 1994]. At lower altitudes, measure- 
ments of N2O and H2O suggest that mixing out of the 
tropics is much more efficient, occurring on a timescale 
of 3 months or less below 22 km [tfinfsa et al , 1994]. 

The above studies emphasized transport from the 
tropics to midlatitudes, but an equally important is- 
sue is the mixing of midlatitude air into the tropics. 
The timescale associated with mixing into the tropics 
is a key factor in assessing the stratospheric effect of 
high-speed civil transport. Model calculations indicate 
that the sensitivity of ozone loss rates to perturbations 
in NO x levels are much greater at 30 km than at 20 
km [Considine et al t 1995]. It is thus important to 
understand the extent to which aircraft exhaust emit- 
ted at midlatitudes may be mixed into the tropics and 
carried to higher altitudes where its effect on ozone is 
significant. Observations of long-lived gases and reac- 
tive compounds NO y and O3 indicate a small, but possi- 
bly important, amount of in-mixing from the extratrop- 
ics [Avallone and Prather , 1996]. Here we examine the 
bulk properties of stratospheric mixing into the trop- 
ics through analysis of long-lived tracer fields obtained 
from balloon, aircraft, and satellite platforms. The ob- 
servations are coupled to a diagnostic formulation for 
tropical tracer budgets which permits estimation of the 
mean rate of isentropic mixing from midlatitudes to the 
tropics in the lower stratosphere. 

Observations 

We focus on two long-lived species, N2O (nitrous ox- 
ide) and CCI3F (CFC-11). Sources for both gases are 
confined to the troposphere, and the sinks are due pri- 
marily to photodissociation in the stratosphere. Global, 
steady state lifetimes are calculated to be 123 and 44 
years for N2O and CCI3F, respectively [Minschwaner et 
al, 1993]. Local photochemical lifetimes at 22 km in the 
tropics are considerably smaller, of the order of 25 years 
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for N 2 0 and 1 year for CC1 3 F in the annual mean. The 
wide range of photochemical time constants between 
these two species is useful for the problem at hand, as 
effects of isentropic mixing are more easily quantified 
if the mixing timescale is shorter than the local pho- 
tochemical lifetime. Another consideration, however, is 
the mean difference between tropical and midlatitude 
mixing ratios along a given isentropic surface. 

Figures 1 and 2 show balloon and aircraft measure- 
ments of tropical CCI3F (Figure I) and N 2 0 (Figure 2) 
spanning 5°S to 10°N by Goldan et ai [1980] and Ved- 
der et ai [1978, 1981]. Midlatitude data shown in Fig- 
ures 1 and 2 are from 41°N by Goldan et ai [1980]. The 
Goldan et al. N 2 0 data are adjusted according to the 
correction factor recommended by Goldan et ai [1981]. 
Uncertainties for these data range from 10% to 25%. AH 
observations were made during the period 1976 to 1979 
and are scaled to 1980 using a height-independent factor 
derived from observed trends in the troposphere [Elkins 
et aij 1993]. The altitude coordinate of the original 
data is mapped to annual mean isentropes in Figures 

1 and 2 using climatological pressure and temperature 
fields [ Anderson et a/., 1986]. 

Curves drawn through the data in Figures 1 and 

2 are least squares, quadratic exponential fits ( y = 
exp (ax 2 + bx -f c)) which are used in the subsequent 
analysis. The limited amount of tropical data precludes 
the use of a nonparametric fit for this work. Choice 
of the quadratic exponential function is guided by the 
expected form of the tropical vertical profile in the ab- 
sence of isentropic mixing, as discussed below. In addi- 
tion, this is one of the lowest order functions that can 
reproduce the correct asymptotic limits (0 at 00, ap- 



Figure 1. CC1 3 F mixing ratios obtained from balloon 
and aircraft measurements by Goldan et ai [1980] and 
Vedder et ai [1978, 1981] mapped to climatological po- 
tential temperature as discussed in the text. Black di- 
amonds are data obtained between 5°S and 10° N; gray 
squares are for data from 41°N. Dashed curves are least 
squares fits to the data. The solid curve labeled “zero 
mixing” is calculated according to (2). 
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Figure 2. Same as Figure 1 but for mixing ratios of 
N 2 0. 


proximately constant mixing ratio below tropopause) 
and has been used previously by Ehhalt et al. [1983] 
to best describe their measurements of tracer vertical 
profiles at midlatitudes. All of our least-squares fits are 
constrained to match tropospheric values [Elkins et al., 
1993] at the tropopause. 

The data shown in Figure 3 consists of CC1 3 F mea- 
surements obtained by the National Oceanic and Atmo- 
spheric Administration (NOAA) Airborne Chromato- 
graph for Atmospheric Trace Species (ACATS-4) on- 
board the NASA ER 2 aircraft [Elkins et al., 1996], The 
observations were made during the Airborne Southern 
Hemisphere Ozone Experiment Measurements for As- 
sessing the Effects of Stratospheric Aircraft (ASHOE/ 
MAES A) (March-November 1994), and the Strato- 
spheric Photochemistry, Aerosol, and Dynamics Expe- 
dition (SPADE) using ACATS-2 (May 1993). Mid- 
latitude data consist of measurements between 30° to 
40°latitude. Discrimination against anomalously high 
or low data influenced by the short-term effects of plan- 
etary wave activity is accomplished using modified po- 
tential vorticity (MPV) [Lait, 1994] as a quasi-conserved 
tracer. The midlatitude data are constrained by values 
of MPV between 8x 10 -6 and 1.6x 10 -5 K m 2 kg -1 s -1 . 
MPV is obtained from daily National Meteorological 
Center analyses and interpolated to the measurement 
location. Tropical CC1 3 F data, defined here for ±15° 
latitude and |MPV| < 6x 10 -6 K m 2 kg~ l s -1 , were ob- 
tained only during ASHOE/MAESA. The uncertainty 
of each measurement is estimated to be < 3%. 

Figure 4 presents N 2 0 measurements obtained from 
the SPADE and ASHOE/MAESA ER 2 flights. During 
SPADE, N 2 0 was measured by the ATLAS tunable 
diode laser instrument [Podolske and Loewenstein, 1993], 
The measurement uncertainty is 5% (lcr). Measure- 
ments of N 2 0 during ASHOE/MAESA are taken from 
the NOAA ACATS instrument. These data have a 
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Figure 3. CCI3F mixing ratios measured during 
ER 2 flights from the Airborne Southern Hemisphere 
Ozone Experiment - Measurements for Assessing the 
Effects of Stratospheric Aircraft (ASHOE/MAESA) 
and the Stratospheric Photochemistry, Aerosol, and 
Dynamics Expedition (SPADE) missions, as a func- 
tion of potential temperature obtained by simultane- 
ous measurements of pressure and temperature. Black 
points are from ±15° latitude with |MPV| < 6 x 10 -6 
K m 2 kg“ l s’ 1 ; gray points are from 30°-40°N with 
8 x 10-* < MPV < 1.6 x 10 “ 5 K m 2 kg " 1 s’ 1 . As 
in Figure 1 , dashed curves are fits to the data and the 
solid curve is the zero mixing solution given by ( 2 ), 

measurement uncertainty of less than 3%. Nearly all 
of the ER 2 data lie below the 500 K isentropic level. 
At higher 9 , the data shown in Figure 4 are from the 
Cryogenic Limb Array Etalon Spectrometer (CLAES) 
[Roche et a/., 1993] onboard the Upper Atmosphere Re- 
search Satellite, measured between February 12 and 
March 16, 1993. Typical uncertainties in these CLAES 
N 2 O measurements are approximately 20%. Discrim- 
ination of all N 2 O data into midlatitude and tropical 
bins is performed on the basis of latitude and MPV, as 
discussed above for CCI 3 F. 

Calculations 

Inputs of CCI 3 F and N 2 O to the stratosphere are 
governed primarily by the cross-tropopause flux in the 
tropics. The local budget of these tracers in the tropi- 
cal lower stratosphere is dominated by effects of upward 
transport, photochemical loss, and mixing with the ex- 
tratropics. Focusing exclusively on the tropics, the con- 
tinuity equation for the tropical mixing ratio x of either 
gas is given by 

d l + i % + L * + ;<*-*'> = ° <» 

where 9 is potential temperature, 9 is the rate of radia- 
tive heating, L is the photochemical loss rate, and x* 



Figure 4. Same as Figure 3, except for N 2 O from 
the ASHOE/MAESA and SPADE ER -2 missions. Data 
above about 500 K are from the CLAES satellite instru- 
ment on board the UARS satellite. 


is the midlatitude mixing ratio. Quasi-horizontal mix- 
ing along surfaces of constant potential temperature is 
parameterized using a mixing timescale r. This is the 
local timescale for mixing into the tropics on a given 
isentropic surface. Our approach assumes that isen- 
tropic mixing between tropics and midlatitudes is slow 
compared with mixing confined solely within each re- 
gion. This point of view is supported by bomb and vol- 
canic debris studies as well as observations of differences 
in tracer correlations between tropical and midlatitude 
regimes [Murphy et a/., 1993]. 

Using the analytic fits to the measured vertical pro- 
files of x> x\ dx/99, and calculated values of 9 and 
L, we can apply (1) to determine r as a function of 9. 
Photochemical loss ( L ) due to ultraviolet photodissoci- 
ation is calculated using the high-resolution radiation 
code described by Minschwaner et al [1993] and in- 
cludes effects of scattering in the 190-210 nm spectral 
region. These loss rates give rise to calculated values 
of global and local lifetimes noted previously. The N 2 O 
loss includes a small contribution (about 10 %) due to 
reaction with 0 ( X D). Long-term atmospheric growth 
rates (dx/dt) are specified from analysis of tropospheric 
measurements (see Elkins et al [1993] for discussion of 
methods). This term is small compared with either ver- 
tical advection or photochemical loss except in the case 
of CCI 3 F data by Goldan et al [1980] and Vedder et al 
[1978, 1981] (hereinafter referred to as Goldan- Vedder), 
where 1980 growth rates were 5 . 7 % per year. 

Radiative heating (0) is calculated using the Commu- 
nity climate Model 2 radiative transfer code [ Briegleb , 
1992a, b] initialized with tropical temperature, ozone, 
and water vapor profiles observed during the Central 
Equatorial Pacific Experiment [Dessler et al t 1996J. 
Tropospheric clouds amounts are specified according 
to annual mean values from the International Satellite 
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Cloud Climatology Project [Rossow and Schiffer , 1991]. 
Longwave fluxes to space and shortwave albedo are in 
good agreement with measurements from the Earth Ra- 
diation Budget Experiment [ Barkstrom , 1984]. Figure 
5 shows the vertical profile of radiative heating. Im- 
plied ascent velocities in the tropical lower stratosphere 
(about 0.025 km day -1 ) are consistent with the mean 
value determined by Holton [1990] on the basis of global 
mass balance. 

Results 

The zero mixing curves shown in Figures 1-4 are solu- 
tions to the homogeneous form of (1) using calculated 
values of 8 and L. That is, we set r = oo and solve 
for x under the assumption of no tropical-midlatitude 
mixing. The vertical profile thus takes the form of an 
exponential 



where 7 = (1 /x)(9x/9t) is the long-term growth rate, 
8 0 is the potential temperature of the tropopause, ap- 
proximately 375 K, and x° is the tropopause mixing ra- 
tio. The argument to the exponential function in (2) for 
both gases is well approximated by a quadratic polyno- 
mial, thus motivating the use of quadratic exponential 
fits to the observed vertical profiles. 

The poor agreement shown in Figures 1-4 between 
the zero mixing solutions and tropical data is a clear 
indication of entrainment of midlatitude air into the 
tropical ascent region. The discrepancies are not related 
to errors in computed photochemical loss: uncertainties 
in photolysis frequencies are no larger than 20% [Min- 
schwaner et al 1993] which cannot explain the differ- 
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Figure 5. Calculated values of 8 used in (1) (solid 
curve, lower axis). Approximate vertical velocities are 
indicated by the dotted curve (upper axis). 


ence between the zero mixing curve and the data. Possi- 
ble errors in computed radiative heating also cannot be 
responsible, since it is easily shown that (2) gives rise to 
a power law relationship between CCI3P and N2O mix- 
ing ratios which is independent of the diabatic heating 
rate: K 

XCCI 3 F — XcCUF ('o 2 ° } (3) 

\an 3 0 / 

where K («22 in the tropical lower stratosphere) is the 
ratio of CCI3F to N2O loss plus growth terms (the nu- 
merator of the integrand in (2) above). We found large 
differences between the power law relation and the trop- 
ical CCI3F-N2O correlations obtained from Figures 1 
and 2, indicating that the zero mixing discrepancies 
cannot be due primarily to errors in computed radia- 
tive heating rates. 

As mentioned above, timescales for isentropic mixing 
between the tropics and midlatitudes can be obtained 
from (1) using a combination of calculated quantities 
and mean vertical tracer profiles (analytic fits) based 
on the observations shown in Figures 1-4. Timescales 
derived in this fashion are shown in Figure 6. Owing 
to the higher measurement accuracies and availability 
of simultaneous observations of potential temperature, 
there is reason to place a higher weight on results from 
the ASHOE, SPADE, and CLAES (hereinafter referred 
to as A-S-C) data. These mixing time profiles are in fact 
bracketed by results from the Goldan-Vedder analyses 
shown in Figure 6. 

Despite large differences in absolute magnitudes, a 
common feature of all data sets is the relatively slow 
mixing between about 430 and 500 K (19-21 km). Faster 
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Figure 6. Timescale for mixing into the tropics based 
on fits to the data shown in Figure 1 (dash-dot, labeled 
G-V (Goldan-Vedder) CCI3F), Figure 2 (dotted, labeled 
G-V N 2 0), Figure 3 (dashed, labeled A-S (ASHOE- 
SPADE) CCI3F), and Figure 4 (solid, labeled A-S-C 
(ASHOE-SPADE-CLAES) N 2 0). Error bars represent 
the level of uncertainty for the A-S-C N 2 0 derived mix- 
ing time, as shown in Figure 7. 
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rates of mixing are indicated above and below this 
implied barrier, although these tendencies should be 
viewed with caution in light of the large level of un- 
certainty (see below). Nevertheless, it is of interest to 
note that the morphology of mixing suggested from Fig- 
ure 6 is similar to that observed for export out of the 
polar vortices [e.g., Schoeberl et a/., 1992]; slow mixing 
in both cases is associated with strong gradients in po- 
tential vorticity [Trepte and Hiichman , 1992; Plumb et 
a/., 1994]. 

Uncertainties in the derived mixing times are esti- 
mated to be about 55% near 450 K, decreasing to 40%- 
46% above 500 K in the case of the A-S-C N 2 0 mixing 
time. Contributions to this uncertainty are displayed in 
Figure 7 which also includes the total uncertainty as- 
suming that individual components add in quadrature. 
Calculated values of radiative heating (±20%-25% in- 
cluding effects of tropical cirrus) and photochemical loss 
(±15%-20% [Minschwaner et a/., 1993]) both make con- 
tributions to the overall uncertainty. Additionally, the 
large scatter of the data introduces uncertainties in our 
\, and dx/98 terms in (1). In regions where the 
quantity \ — \ f is comparable to the uncertainties in \ 
and \ 7 , such as below 400 K, computed mixing times 
are uncertain by at least a factor of 2. 

Owing to the large temporal range of the data, our 
results represent mean values averaged over at least 1 
year. These mixing times thus represent the bulk ef- 
fect of entrainment into the tropical lower stratosphere 



UNCERTAINTY IN MIXING TIME SCALE (%) 

Figure 7. Components to the uncertainty for the A-S- 
C N 2 0 mixing timescale shown in Figure 6. The curves 
labeled 6t\ and Sr 2 are contributions from uncertainties 
in \ and respectively, assuming these uncertainties 
both increase with 9 from 7 to 15 ppb within the alti- 
tudes of interest; St 3 is due to an 18% uncertainty in L\ 
8t 4 arises from a 50% uncertainty in dx/d6 at 380 K, 
decreasing to 20% at 650 K; and Sr 5 originates from a 
23% uncertainty in 9. Total uncertainty is calculated by 
adding the individual components in quadrature. Un- 
certainties for the other mixing profiles shown in Figure 
6 are generally larger. 


due to synoptic and planetary scale wave mixing. The 
strength of the subtropical barrier may also depend on 
season as well as the phase of the quasibiennial oscil- 
lation in tropical winds and temperature [ Hiichman et 
a/., 1994], but these effects are beyond the scope of this 
paper. Model simulations also suggest hemispheric dif- 
ferences in mixing out of the tropics [ Chen et a/., 1994]. 
We found small hemispheric differences in the CLAES 
midlatitude distributions, although the impact on com- 
puted mixing into the tropics was not significant. 

Species used in this analysis must satisfy two crite- 
ria. First, the local photochemical lifetime in the tropics 
must be larger than the derived mixing time so that the 
effect of mixing is apparent in the tropical profile. Sec- 
ond, there must be a clear distinction between tropical 
and midlatitude profiles (given by x X* * n (!))♦ We 
believe N 2 0 is better suited for this analysis because (1) 
the local photochemical lifetime exceeds 1 year in the 
tropics below about 28 km and (2) the reaction with 
O l D becomes more important at middle and high lat- 
itudes in setting the total loss rate, thus producing a 
sufficiently large gradient between tropical and midlat- 
itude profiles. Our highest confidence therefore lies in 
results from the A-S-C N 2 0 data. It is not clear why 
timescales derived from the Goldan-Vedder N 2 0 data 
are significantly larger above 500 K. 

A quantity of interest is the fraction of midlatitude 
air entrained into the tropical ascent region. If (I) 
is rewritten in terms of the difference A between the 
tropical and midlatitude mixing ratio of an infinitely 
long-lived tracer in steady state ( L = 0, d/dt = 0) and 
we define £ = (A c — A)/A 0l where A 0 is the initial dif- 
ference at the isentropic surface intersecting the tropical 
tropopause, then £ is determined by 

<(«) = l-exp (y- J 1 ®) 

Owing to the altitude dependence of mixing and tracer 
vertical profiles, this quantity is not the same as the 
apparent mixing fraction derived along surfaces of con- 
stant 9. For example, the apparent mixing fraction (ob- 
tained by simply comparing the observed and zero mix- 
ing profiles on a constant theta level) for the A-S-C N 2 0 
data is about 30% at 650 K. As shown below, the ac- 
tual fraction of midlatitude air is approximately 65%. 
The reason this value is much higher than the apparent 
mixing fraction is that much of the mixing is occur- 
ring at lower altitudes where the difference between the 
tropical and midlatitude mixing ratios is smaller. 

Figure 8 shows the calculated profiles of £ correspond- 
ing to mixing timescales from Figure 6. Differences be- 
tween data sets are less than 40% above 500 K despite 
much larger disparities in profiles of r; apparently, £ is 
a more robust quantity for the present method of analy- 
sis. As an indication of the effect of uncertainties in the 
vertical profile of r, Figure 8 also includes a £ profile 
corresponding to a fixed r of 12 months. The fraction of 
midlatitude air entrained into the tropics increases from 
0% at the tropopause to more than 50% above 550 K 
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Figure 8. Fraction of midlatitude air entrained into 
the tropics based on the mixing time curves shown in 
Figure 6 (see (4) and discussion in text). Additionally, 
a C profile is shown that corresponds to a fixed r of 
12 months, independent of 9. Error bars represent the 
level of uncertainty for the mixing fraction derived from 
the A-S-C N2O analysis which includes the uncertainty 
in r and in 9. The latter uncertainty, however, very 
nearly cancels in the expression for (. 

for all cases except the ASHOE-SPADE CCI3F profile, 
which reaches 32% at 490 K. Note that the definition of 
midlatitude used here does not distinguish between the 
ultimate origins of midlatitude air parcels, only that 
they have spent sufficient time in the extratropics to 
acquire midlatitude characteristics of long-lived tracer 
amounts. 

Our results for ascent with quasi-continuous mixing 
are qualitatively consistent with observations of the up- 
ward propagating seasonal signal of 2xCH4 + H2O pre- 
sented by Mote et al. [1996]. They found that the 
amplitude of the seasonal signal (comparable to 1 — £ 
in the present analysis) diminished sharply between the 
tropopause and about the 500 K level but changed much 
less rapidly at higher altitudes, in agreement with the 
mixing fractions shown here in Figure 8. The signal 
attenuation below 500 K implies a lower limit to the 
mixing timescale of 7 to 9 months, also in agreement 
with our values of r shown in Figure 6. However, the 
lower limit of 15-18 months derived by Mote et al. be- 
tween 500 and 700 K is somewhat longer than our best 
estimate given by the A-S-C N2O data in Figure 6. 

Conclusion 

The distributions of long-lived tracers N2O and CCI3F 
in the lower stratosphere indicate that a potential bar- 
rier exists for isentropic mixing of midlatitude air into 
the tropics. The tropical stratosphere may be effec- 
tively isolated for more than 1 year between 430 and 
500 K. Although enhanced rates of mixing are inferred 


above and below this region, the level of uncertainty 
does not preclude the possibility of a fixed value of 
r, independent of 9. Additional data from the tropi- 
cal lower stratosphere, for example, from the ongoing 
Stratospheric Tracers of Atmospheric Transport mis- 
sions, would be highly desirable to more effectively con- 
strain mean tropical profiles of long-lived tracers. It is 
unlikely, however, that the present method can be used 
to unambiguously derive a vertical profile of mixing in 
which variations with respect to 9 will exceed the level of 
uncertainty. The mean timescales derived here suggest 
that further analysis would be best suited using N2O or 
molecules with similar photochemical time constants, 
such as CCI2F2 (CFC-12). Recent work by Volk et al. 
(In situ measurements constraining exchange between 
the tropics and the middle latitudes in the lower strato- 
sphere, submitted to Science, 1996) has provided mean 
mixing times below 490 K by employing tracer corre- 
lations and involving a more extensive suite of source 
gases and photochemically active species. Their mean 
in-mixing time of 13.6 months below 490 K agrees fa- 
vorably with values reported here. 

The implied entrainment of midlatitude air into the 
tropical ascent region is remarkably large; up to 50% 
of the air above the 550 K surface is of midlatitude 
origin. Recently, Plumb [1996] has described a concep- 
tual “pipe” model of the tropical ascent region which 
successfully addresses shortcomings in the Plumb and 
Ko [1992] global mixing formulation. We find evidence 
here which supports a pipe-like behavior in the lower 
stratosphere. Owing to the integrated effect of isen- 
tropic mixing, however, the tropical import barrier ap- 
pears to possess an altitude-dependent permeability. 
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